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Please enter short abstract We are exploring the possibility of extending the cycle life of 

battery systems by using a charging profile that minimizes cell 
damage. Only nickel-hydrogen cells are discussed at this time, but 
applications to lithium-ion cells are being considered. The process 
first requires the development of a fractional calculus based 
nonlinear dynamic model of the specific cells being used. The 
parameters of this model are determined from the cell transient 
responses. To extend cell cycle life, an instantaneous damage rate 
model is developed. The model is based on cycle life data and is 
highly dependent on cell voltage. Once both the cell dynamic model 
and the instantaneous damage rate model have been determined, the 
charging profile for a specific cell is determined by numerical 
optimization. Results concerning the percentage life extension for 
different charging strategies are presented. The overall procedure is 
readily adaptable to real-time implementations where the charging 
profile can maintain its minimum damage nature as the specific cell 
ages. 


This is a preprint or reprint of a paper intended for presentation at a 
conference. Because changes may be made before formal 
publication, this is made available with the understanding that it will 
not be cited or reproduced without the permission of the author. 
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Control Philosophy 

Two Phases of Control System Design : 

Phase I is the design of optimal trajectories and associated inputs, that 
move a given plant from one operating condition to another, while 
minimizing some performance measure. Requires a nonlinear dynamic 
model of the specific system. 

Phase II is the design of a trajectory following controller (sometimes 
called a regulator or tracker) that provides a real-time control input 
perturbation to keep the plant operating near the designed optimal 
trajectory. Usually uses a linearized dynamic model of the specific system. 

Phase I: 

Our control philosophy is to charge the NiH2 cell in such a way that the 
damage incurred during the charging period is minimized, thus extending 
its cycle life. Requires nonlinear dynamic model of NiH2 cell and a damage 
rate model. We must do this first. 

This control philosophy is generally considered damage mitigating control 
or life-extending control. 
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Overview of NiH2 Electrochemistry 


+ 



Normal charge-discharge operation of a nickel-hydrogen cell: 


positive nickel electrode: NiOOH + H 2 0 + e ' Ni(OH) 2 + OH ' 

negative hydrogen electrode: j H 2 + OH ' t ^ * H 2 0 + e ' . 
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Essentialized performance model of NiH2 cell. 


storage 


diffusion 



Jiicht irgi/ 

S + D + e " ' reaction products 
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Electrode Behavior: Faraday’s Law 

• For the discussion of the next two sections, consider the electrode process : 

aA + bB < > cC + dD + ne~ m 

• This equation represents both a chemical process and an electrical process. 

• The reaction rates can be completely determined using the electrical process by seeing 
that the chemical conversion can only occur if electrons are either arriving or leaving. 

• Thus, the chemical conversion rates are controlled by, or measured by, the electrical 
current passing through a given electrode. 

• Then recognizing that the rate of electron production is related to electrical current /, 
the following rate equations result : 

•= C !l£ < L c L 0 - c dC _ | a dA _ ' b dB 

dt n dt n dt n dt n dt 

where <I e - is the charge of a single electron. 
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Electrode Behavior: Electrode Equation 

Consider the fluxes, y, of the species in the forward and reverse reactions 
at the electrode, 

* = j f ~ Jr - 

Assuming that the species fluxes are proportional to concentrations, yields 

i = k f c A c B — k r c c c Da 

The rate constants, can be related to the electrical potential across the 
electrode-electrolyte interface using free energy considerations, 

k, = j 

k = ^ oe (F//?r){-a){v-v fl ) 


Inserting these gives the electrode equation. 


: = *0 { C A C B ' 




c d (F i RT)(-a)( 
L C c D * 


j 


The approach used is often referred to as the Butler-Volmer approach. 
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Linear-in-the-parameters Electrode Equation 

We propose a linear-in-the-parameters approximate solution to the electrode equation; 


' = k { + k 2 In (l +| / |)sgn(/) + k 3 In (c d )+ k A ln(i - c s ) 


where the fc’s are parameters to be determined from data. 

This represents a compromise between the Tafel and the Nernst solutions of 
the electrode equation: 

The Nernst solution assumes that the current is so small as to be negligible 


RT 

nF 


In 




\ c d ) 


The Tafel solution assumes that the current is large in one direction or 
the other, which means that one of the two exponential terms is negligible 

v= v n + — -\n{k 0 ) -ln(i ) m 


anF 


anF 
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Essentialized Model Overview 


Terminal behavior: 


current into the battery is +/, 

the terminal voltage is +v, 


stored material with self-discharge: 


dc(t) v x 1 . . 

-^-’ = ,(0-— c,«* 


\%d 


diffusing material: 

C d (0 ~ C bulk — A q K t ) , 

electrode-equation: 


v(c, , Cj , /) = k l +k 2 ln(l + 1 / |)sgn(/) + k 3 In (c d ) + k 4 ln(l - c s ) . 
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Parameter Determination 


The cell chosen is the NSWC Crane Pack ID 3602G (Gates): 
rated at 65 AHr 
uses 3 1 % KOH concentration 
maintained at 1 0 degrees C 
charge-discharge profile is a square wave 
with 35% depth-of-discharge (DOD) 

104% recharge ratio 

current is 26.29 A for 54 minutes charging 
-37.92 A for 36 minutes during discharge 
note that 65 AHr - 3900 AMin, 35% of 65 AHr = 1365.1 
AM in . 
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Essentialized Model with Identified Parameters 


Terminal behavior: current into the battery is +i, 

the terminal voltage is +v, 

stored material with self-discharge: 

= i(t) - 0.0002085c 5 (?) 
dt 


diffusing material: 

c d (t)= 1 -0.001036 0 ^' 903 7(0, 

eiectrode-equation: 

v =1 ,3656+0.02651n(l + | / |)sgn(/')+0.02291n(c^) 

-0.0262 ln((l .005 * 3 900 - c s )/ 3900) 
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Damage Mechanisms for NiH2 Batteries 


Formation of y- phase NiOOH: 


overcharge 

P -NiOOH ► y -NiOOH 


i 

i 


i 

V 

charge 


discharge 

charge 




f 


▼ 


discharge 


p -Ni(OH) 2 < a -Ni(OH) 2 

standing in KOH 

Bode’s Solid phase relationships for a NiOOH electrode. 
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Damage Mechanisms for NiH2 Batteries 

Formation of 02: 

Overcharge: 1 ) continuing to charge the cell after all the p -Ni(0H)2 
has been converted 
2) the charging current is too large 

The effect of this is the formation of 0 2 at the nickel electrode, 
along with heating. 
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Damage Mechanisms for NiH2 Batteries 


Heating: 

Sources: 1) heat of reaction 

2) formation of 0 2 

3) electrical current 

Results in: 1) the formation of r-NiOOH: 

which a) reduces cell capacity 

b) does physical damage to the cell 

Results in an increase in self-discharge reaction rates. 
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Battery Continuum Damage Modeling 

• Many Possible Damage Mechanisms 
•Hard to model all these 


• Overall Birth to Death Data will be Used Instead 


• Crane Database Provides much Information 


• Green-Hoffman Data Taken as Starting Point 
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Cycles to Failure 


Battery Continuum Damage Modeling 


Green-Hoffman DataNiH 2 atT = 10 °C 

DOD C ^ cles 

To 

Failure 
35% 38,000 

50% 19,000 
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Green-Hoffman Battery Life Model 
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Continuum Damage Model 


Based on G-H Data 

N f =1885.04/- 621(l-DOZ)) 


= 19151 1.73e“ 4621DOD , at 1 0° C 
For constant damage per cycle 


D = 


1 


eye 


N 


= 5.222 x\0~°e 


6 4.621 DOD 


fcH 

DOD = c,v 0 


D cvc = 5.222 xlO -6 e 4 ' 62lci, ’“ 
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Continuum Damage Mode! 


jS(v)dv = D cyc = 5.222x\0~ 6 e 4 ' 621c '‘“ 

cycle 


where S(v) = C ^~ 

dv 

For damage on charging only 


voltage referred 

damage rate 


* max 

p(v)flfv = =5.222xl0 _6 e 4 ' 621C| "* 

* min 
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Continuum Damage Model 


v —v 
max min 

p(y + v min ) dv = 5.222x 1 0” 6 e 4 62 lc > v " 

0 

= 5.222 x 1 0 -6 e 4 ' 621c ‘ } 

Thus inferring 


<?(v + v m J £ 5.222x10^4.621 c, e 4 ' 621c ' 1 ’ 
= 2.413 lx I0~ 5 c 1 e 462lc ''' 

Hence 

<y(v) =2.4131xl0“ 5 c, e 4 - 621c ' (v “‘“ ) 
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Continuum Damage Model 

Instantaneous Damage Rate 

• , . dD dD dv x. ...... 

d(<) ~~J7 = ~r^r" S(vU))v{t) 

dt dv dt 


Requiring Positive Damage 


23(0 = — = — 

dt dv 


23(0 = 2.4131x10' 


dv 

dt 


c x e 


= #X0) |v(/)| 

4.621 c, (v— v mln ) |y(/)| 


23(0 = 2.4131x10 _5 c, e 4 ' 621 c ><'"- , - 2 > |v(,)| 
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Modified Continuum Damage Model 


For zero damage at zero POD 

D cy c=c 2 {e c,D0 D -e c ' n ) 

l)=J- 

c 2 and c 3 are determined to match G-H Data 

D = 1.0404xl0“ 5 b 3M2DOD -l)= — 

' V / 

J 

Repeating the previous process gives 

D(t) = 3.7475x10" 4 c 1 e 3Mlc ^~'- 2 '> |v(;)j 



Battery Life Models 
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Modified Continuum Damage Model 


In terms of current 


No- 


dD _ dD 

dv 

dq 

= S(y(t)) 

dv 

dt dv 

dq 

dt 


dq 



where — is slope of charging curve 

dq 

For zero damage when voltage rate goes negative 

D{t) = 3.7475xl0“ 4 c 1 c 36O2qM - 2) v(O, v(/)>0 
D(t) = 0, v(/)<0 
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Control Philosophy 

Two Phases of Control System Design : 

Phase I is the design of optimal trajectories and associated inputs, that 
move a given plant from one operating condition to another, while 
minimizing some performance measure. Requires a nonlinear dynamic 
model of the specific system. 

Phase II is the design of a trajectory following controller (sometimes 
called a regulator or tracker) that provides a real-time control input 
perturbation to keep the plant operating near the designed optimal 
trajectory. Usually uses a linearized dynamic model of the specific system. 



Phase I: 

Our control philosophy is to charge the NiH2 cell in such a way that the 
damage incurred during the charging period is minimized, thus extending 
its cycle life. Requires nonlinear dynamic model of NiH2 cell and a damage 
rate model. Now that we have this we can begin the control design process. 
The specific control philosophy is employed is generally considered damage 
mitigating control or life-extending control. 
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Performance Measure 


The performance measure to be minimized is the accumulated 
damage per recharge cycle: 

j = w fs k {t f )~ c s (t f )f + dt 

0 

c](t f ) is the desired stored charge at the end of charge, 

w fs is the cost weighting , 
dDjt) 

j t is obtained from the damage model, 
tf for our problem is 54 Min. 
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65 %- 100% recharging, voltage profile 



Comparison of 65%-100% charging methods, voltage 
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65%-100% recharging, damage profile 



Comparison of 65%-100% charging methods, damage. 
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65%-100% recharging, current profile 
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Comparison of 65%-100% charging methods, current. 
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60%-95% recharging, voltage profile 



Comparison of 60%-95% and standard charging methods, 
voltage. 
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60%-95% recharging, damage profile 



Comparison of 60%-95% and standard charging methods, 

damage. 
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60%-95% recharging, current profile 
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Optimal Charging Summary 



n 

Cycles to 
Failure 

% Life Extension 
relative to constant 
current 

% Life Extension 
relative to constant 
current-plus-taper 


0.000039269 

25465 

0% 

not applicable 


0.000026316 

38000 

49.22% 

0% 

Life Extending 
Charging 
65%-100% Cycle 
of Figure 5.4 with 
abs(dv/dt) damage 
rate 

0.000020780 

48123 

88.98% 

26.64% 

Life Extending 
Charging 
60% -9 5% Cycle 
of Figure 5.5 with 
abs(dv/dt) damage 
rate 

0.000019535 

51190 

101.02% 

34.71% 

Life Extending 
Charging 
65%-l 00% Cycle 
of Figure 5.6 with 
only +dv/dt damage 
rate 

0.000022080 

45290 

77.85% 

19.18% 


Comparison of damage for various charging methods, assumes damage only during charge. 
Percentage life extension is reduced proportionately for damage during discharge. 


Electrical and 
Computer 
Engineering 































Control Design - Phase II: Tracking 


Input Measured output 



Real-time observer structure. 
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Control Design - Phase II: Tracking 












Summary 


• Control Philosophy 

• Essentialized Model Development 

• Damage Model 

• Optimal Life-Extending Charging 

• Tracking Controller 

• Real-time Parameter Identification Development 

• Application to Lithium based cells 
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